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Abstract 

We consider a model of a molecular junction made of BDT (benzene dithiol) molecule trapped between 
two Au(100) leads. Using the ab initio approach implemented in the SIESTA package we look for the 
optimal configuration of the molecule as a function of a distance between the leads. We find that for the 
distance long enough the energy of the system is minimized when the molecule is bonded asymmetricaly, 
i.e. chemisorbed to one of the leads, whereas for the distance shorter than 12 A the energy is minimized 
for the molecule sitting in the middle between the leads. We discuss possible consequences of the above 
findings for the transport properties of the junction. 

PACS numbers: 73.63. -b, 85.65. +h, 73.40. -c 
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Molecular junctions made of two leads bridged by a single molecule are intensly stud- 
ied recently both experimentally and theoretically as basic units in the systems of molecular 
electronics^^. Unlike the traditional silicon based devices, the electronic transport through the 
molecular devices very much depend on subtle details of molecular configuration between the 
leads, due to importance of interference effects at the nanoscale. In several ab initio works the 
influence of a type of bonding between the metal surface and the bridging atom of the molecule 
(on-top, hollow place, bridge contact structures) on the current-voltage (I-V) characteristics was 
thoroughly studied for both Au(lll) and Au(100) surfaces 4 -^. Moreover the variability of the 
transport properties of the junctions with the change of width of the leads (being also the kind of 
interference phenomenon) was demonstrated 4 . So far however a full understanding of the relation 
between the distance between the leads and the details of the I-V characteristics seems to be absent 
in the literature. 

In this paper we study the evolution of electronic and transport properties of a molecular junc- 
tion with change of the distance between the leads. We consider benzene dithiol (C6H4S2 or BDT) 
molecule placed between the gold leads, with the sulfur atoms bound to the Au(100) surface at 
the hollow positions, i.e. at the centers of squares formed by surface Au atoms. We first look for 
the optimal configuration of the molecule between the leads and analyze how this configuration 
depends on the distance, and next we compute the transmission function and the I-V dependence 
for the optimized geometry of the junction. 

To find the optimal configuration of the Au(100)-BDT-Au(100) junction we apply the density 
functional approach as implemented in the SIESTA package 7 . In our SIESTA computation we 
consider a periodic 3D lattice with a unit cell consisting of BDT molecule and fragments of the 
opposite leads, which include total of 144 Au atoms (see Figure[l]). 

Such a lattice represents well enough the typical experimental setup as used in the mechanically 
controlled break junction (MCBJ) experiment 1 which involves a single BDT molecule provided 
a number of molecular orbitals of BDT is small enough as compared to a number of Au atoms. 
Moreover, we have to avoid overlapping between the basis functions of BDT molecules from the 
neighbouring unit cells what means that the transverse size of the unit cell should be greater than 
the sum of the transverse size of the molecule and the maximal diameter of the basis function. 

All the numerical results presented below are obtained using local density approximation 
(LDA) for the exchange-correlation functional of the DFT method 9 . The results of other works 
as well as our own calculations show that neither nonlocal correction (i.e. GGA) nor a spin- 
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FIG. 1: A perspective view of the optimized unit cell of the Au(100)-BDT-Au(100) molecular junction 
obtained from gOpenMol for di ea d=ll A. 

-dependent functional (LSDA) change the results substantially. We use the Troullier-Martins 8 
pseudopotential to represent potential of the atomic cores. Our function basis is restricted to 
single-zeta (SZ) functions and we checked some results using the DZP basis. The radii of all 
the pseudoorbitals were determined by the pseudo-atomic orbital energy shift equal to 0.02 Ry. 
With the above restrictions all our electronic structure and transport computations (except the ones 
using the DZP basis) were possible using 2GB of a computer memory. We find the optimal con- 
figuration of the junction with the conjugated gradient method, relaxing also the positions of 4 Au 
atoms binding the molecule in each lead. A locally stable configuration is found numerically if the 
forces acting on each of the relaxed atoms are less than 0.1 eV/A. 

For a small distance between the leads, di ea d ~10 A, we find that the molecule centered at 
the middle between the leads attains the local energy minimum. Computations with the molecule 
shifted uniformly out of the center lead to stable configurations with higher energy (or the molecule 
shifted back to the center by the minimization procedure). The energy increases steadily with 
displacement from the central position. We thus conclude that for a small distance between the 
leads the central position represents the global minimum and there is no any other energetically 
stable configuration. The increase of di ea d reduces the overlap of the S orbitals of BDT and Au 
orbitals from the leads surface and the energy of the system increases (see Figure The central 
position remains the global energy minimum until a critical distance di ead ~ 12.25 A, where we 
find that the energy of the system can be reduced by shifting the molecule towards either one of 
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FIG. 2: (a) Dependence of the energy of the junction as a function of the distance between the Au leads. 
Solid line: BDT molecule at the central position, dotted line: BDT molecule shifted to one of the leads. For 
the shifted molecule the distance between the S atom and the lead surface was about 1.75 A. £ =-13100 e V 
is a reference energy value, (b) Energy of the stable configuration of BDT molecule as a function of distance 
between the sulfur atom and the nearest Au atoms for di ea d =12.5 A. The solution for ds-Au ~3.3 A corre- 
sponds to the central position of the molecule. In the both figures circles show locally stable configurations. 



For the lead distance bigger than the critical one we find that the global energy minimum is 
obtained for the binding sulfur atom at about 1.75 A from the nearest lead surface, corresponding 
to the nearest S-Au distance about 2.7 A. The energy difference between the central position and 
the displaced one rises with increasing the distance between the leads and for a very large di ead 
reaches the value ~2 eV. 

For the optimized junction configurations we computed transmission and I-V characteristics 
using the SMEAGOL package 10 which is an overlay on the SIESTA to allow computations of the 
transport in nanostructures. The details of the applied method were described in a recent paper by 
Rocha and coworkers. The current J was computed using the standard equation of the Landauer 
theory^: 

J = J dET(E, V) [f L (E) - f R {E)\ (1) 

where T(E) denotes the transmission function and fa, fu are the Fermi functions corresponding to 
the left and the right leads. The transmission function is computed using non-equilibrium Green 
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function method with the quasiparticle LDA Hamiltonian. For the small distance between the 
leads (<i; ea rf=10.75 A) energy dependence of T(E) exhibits a pseudogap region in the vicinity of 
the equilibrium Fermi level and initially the current shows a relatively slow linear increase with 
the voltage (see Figure 13). With increase of the voltage the source-drain voltage window extends 
first over the region of the HOMO level of the BDT, and next also the LUMO level, giving rise to 
a steer 




E [eV] voltage [V] 

FIG. 3: Transport properties of the junction with di ea d=l0.75 A. (a) Transmission as the function of energy 
for V=0 (thick black curve) and V=2Y (thin curve, red in online edition). The broken lines show the 
corresponding voltage windows, with E = being the position of the equilibrium Fermi level, (b) Current 
as the function of the voltage (line) and the voltage derivative of the current (circles). 

The overall behaviour of the transmission and the current are quite similar to the results of 
the papers of Xue and Ratner 11 , as well as the Ke et al. 4 . As in the cited papers, the current 
at V=2V is greater than the experimental value by at least order of magnitude which may be 
due to too poor treatment of electron correlations by the LDA method. On the other hand, the 
equilibrium conductance obtained here is several times smaller than the one computed in some 
other recent papers(see e.g. 5 ) and also the shape of our differential conductance is more similar to 
the experimental one 1 . We speculate that the difference may be due to the greater transverse size 
of the Au leads in our paper. 

In Figure |U we show the equlibrium transmission functions in the lowest energy configura- 
tions for several values of the distance between the leads. Until di ea d ~12.5 Athe transmission 
changes gradually with the distance. For a small di ea( i the maximum of T(E) at the HOMO peak 
is smaller than unity (T ~0.9) and separated from E F by about ~1.5 eV. With increase of di ea d, 
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the peak shifts towards E F , and finally it splits into two peaks, where T(E) is very close to unity. 
This behaviour resembles the corresponding one of a simple tight-binding model of a two-atom 
molecule. There, for a coupling between the molecule and the leads stronger than some critical 
value (depending on an interatom hopping), the two peak structure merges into a single peak, 
which decreases with further increase of the coupling 12 . This suggests that the transmission of 
BDT can be partly understood in terms of the two-atom model, in which the internal part of BDT 
(all the atoms except the sulfur ones) can be approximately represented by a small and weakly 
energy dependent hopping between the external atoms binding the molecule to the leads. 



FIG. 4: Transmission as the function of energy for V=0 for variable distance between the leads: d; ea d=10 A, 
10.75 A, 11.5 A, and 12 A (the line thickness decreases with the distance). The broken line near E = 
shows the transmission for e^ ea d=12.5 A, with the BDT molecule attached asymmetrically to one of the 
leads in the lowest energy configuration from the Figure|3b). 

As a result of the upward shift of the HOMO peak, the transmission at E F steadily increases 
with the distance between the leads from T(E F ) ~ 0.05 at e^ earf =10 A to T(E F ) ~ 0.15 at 
di ea d=l2 A. When the distance crosses the critical value the central position of BDT molecule 
between the leads is no longer energetically favoured. The transmission in the asymmetric position 
is considerably smaller (see the broken curve in Figure |4j), since it is limited by a very weak 
coupling with the more distant lead. As a result, the current in the asymmetric position is smaller 
than ~0.2 fiA for any voltage. 

In relation to an experiment the last results suggest, that with a carefully controlled increase 
of the distance between the leads one would observe first a steady rise of the transmission. At 
a distance slightly bigger than 12.25 A the current measured at some constant voltage would 
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suddenly drop. On the other hand, a large value of current observed for a bigger distance would 
suggest that a significant reorganisation of the junction takes place. In fact on the basis of the 
molecular dynamical simulation of Kriiger et alM pulling of a string of Au atoms by a molecule 
chemisorbed to an Au(l 1 1) lead can be expected. 

Finally the general comment is in order concerning the value of the computed current, which 
greatly overestimates the ones observed experimentally. So far this seems to be the common 
deficiency of all the DFT and Hartree-Fock 14 based approaches, at least for a weak coupling or 
relatively small molecules. It remains to be seen if a full selfconsistent implementation of a better 
treatment of the electron correlation (e.g. SIC 15 ) will bring a substantial progress here. Our results 
indicate also, that a better understanding of the role of electron correlations in the transport through 
a simple two-atom molecule can be very helpful in this respect. 
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